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ABSTRACT 



Context. The majority of cataclysmic variables observed in the hard X-ray energy band are intermediate polars where the magnetic 
field is strong enough to channel the accreting matter to the magnetic poles of the white dwarf. A shock above the stellar surface heats 
the gas to fairly high temperatures (10-100 keV). The post-shock region cools mostly via optically thin bremsstrahlung. 
Aims. We investigate the influence of Compton scattering on the structure and the emergent spectrum of the post-shock region. We 
also study the effect it has on the mass of the white dwarfs obtained from fitting the observed X-ray spectrum of intermediate polars. 
Methods. We construct the model of the post-shock region taking Compton scattering into account. The radiation transfer equation 
is solved in the plane-parallel approximation. The feedback of Compton scattering on the structure of the post-shock region is also 
accounted for. A set of the post-shock region model spectra for various white dwarf masses is calculated. 

Results. We find that Compton scattering does not change the emergent spectra significantly for low accretion rates or low white 
dwarf masses. However, it becomes important at high accretion rates and high white dwarf masses. The time-averaged, broad-band 
X-ray spectrum of intermediate polar V709 Cas obtained by the RXTE and INTEGRAL observatories is fitted using the set of computed 
spectral models. We obtained the white dwarf mass of 0.91 ± 0.02 Mq and 0.88 ± 0.02 Mq using models with Compton scattering 
taken into account and without it, respectively. 
Conclusions. 

Key words, radiative transfer - scattering ~ methods: numerical - (stars:) white dwarfs - stars: atmospheres - X-rays: stars 



1. Introduction 



Cataclysmic variables (CVs) ar e close binary systems with a 
white dwarf (WD) as a primary (IWarneijri995h . A WD accretes 
matter from a companion star (typically a red dwarf) that fills 
its Roche lobe. This matter forms an accretion disc inside the 
Roche lobe of the WD. In intermediate polars (IPs), the signifi- 
cant magnetic field (10^-10^ G) can disrupt the disc at some dis- 
tance from the WD surface. The matter is then free-falling along 
the magnetic field lines onto the WD and forms a strong shock 
above its surface (.AizuLl973 ). The temperature in the post-shock 
region (PSR) can be very high (~10-100 keV) and the plasma 
there is cooled mainly vi a optically th in bremsstrahlung radi- 
ation in the X-ray band (Fabian et al.l [1976; Lamb & Masteri 
[197^; King & Lasota|ll979r Cooling via Compton scattering is 
less significant due to a relatively low radiation energy density in 
PSR. This contrasts with the accretion columns of X-ray pulsars, 
where Comptonization is the dominant cooling mechanism (see 
[Becker & Wolff 2007, for more details). 

The theo ry of the PSR has been d e veloped in a number 
ofjpapers (Aizul 119731; IWuetal.1 Il994t IWoelk & BeuermannI 



magnetic funneling (not the standard plane-parallel approxima 
tion). The plasma temperature in the PSR depends on the free- 
fall velocity at the WD surface; therefore, the X-ray spe ctra of 
IPs can be used to determine the WD mass (Rothschild et al 



1981b. Such studies have b een perfo rmed bv ICropper et al 
( 1998'). 'Cropper et alj il99^, iRam sav' ("2000^. 'Beardmore et al. 
( 2000 )^ Revnivtse v et all 0004 ). Suleimanov et al. (2005^ 
and [pala nga et al. (2005), using various X-ray data sets. 
ISuleimanov et al.. (2005) demonstrated that, to accurately mea- 
sure the maximum post-shock temperature to infer the WD mass, 
the spectra at high energies (up to at least 100 keV) are needed. 

The IPs are the most luminous and the hardest X-ray 
sources among accreting WDs. The interest in this class of 
objects has been growing in the past few years as IPs have 
been recently proposed as the dominant X-ray source popu- 
lation detectednea£jheGal^cti£_ce^^ Chandra observa- 
tory dMuno et all 120041; iRuiter et alj l2006l) . The IPs also con- 
tribu te significantly to the X-ray diffuse Galactic ridge emis- 
sion dRev nivtsev et al.' '2006"), as well as dominate the hard X-ray 
sky there ( Krivonos et al. 2007a). Moreover, most of the known 



1996 ; Cropper et al. 1999[). The mos t recent inve stigations by 



Canaile et all (i2005[) and lSaxton et~al] (l2005l l2007h have studied 



CVs d etected by INTEGRAL (Bird et al.|[2007[; iKrivonos et al 
l2007b|) and Swift satellites are IPs (see e.g. iBarlow et al 



the role of the two-temperature plasma and considered the dipole 



_ _ --- - - - - - , - ■ 

200l TBonnet-Bidaud et alj|2007t Ide Martino et al.ll2008h . Inflie 
strongly magnetized {B > 10^ G) polar systems, cyclotron radi- 
ation is an important cooling mechanism, which suppresses the 
high temperature bremsstrahlung emission, whilst it should be 
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where Aj^(T) is the cooling function, /ig the electron num- 
ber density, and rii the ion number density. We take Aj^iT) 
for solar chemical composition as calculated and tabulated by 
[Sutherland & Dopita ( 1993). 

The term 1 + A is the amplification factor due to 
Comptonization. In contrast to the usually used integral am- 
plification factor (which is just the ratio of the emergent 
flux to the incident flux of some seed soft photons, see 
llUarionov & Sunyaevll975l:ISvenssonll984l) . we use here the lo- 
cal amplification factor, because we are interested in the verti- 
cal structure of the accretion column. In a homogeneous isother- 
mal slab, the integral and local factors are nearly equal. As the 
Compton cooling rate depends on the local plasma tempera- 
ture and density in a different way compared to the optically 
thin, thermal radiation, the local amplification factor (1 -H A) is 
a function of the geometrical coordinate z- The cooling due to 
cyclotron radiation is not considered here. 

Equation ([1]) has the integral 



Fig. 1. A schematic view of the accretion flow onto the WD 
magnetic polar cap. The accreting gas falls quasi-radially to- 
wards the magnetic poles, forming a standing shock, below 
which the X-rays are emitted. 

negligible for the IPs. This could explain why the majority of 
the CVs observed in the hard X-ray band are IPs. 

To better understand the hard X-ray production in IPs, one 
needs to model their observed spectral properties. The goal of the 
present paper is to investigate the influence of Compton scatter- 
ing on the structure of the PSR and the emergent spectrum. We 
developed the physical model of the PSR and compare the com- 
puted spectral models to the broad-band X-ray spectrum of the 
intermediate polar V709 Cas obtained by RXTE and INTEGRAL 
observatories. We have also determined the WD mass in this sys- 
tem. 

2. Model of the post-shock region 

2.1. Main equations 

The structure of the stationary post-shock region (the accretion 
and emission geometry is illustrated schematically in Fig.[T]) in a 
plan e-parallel, one-dime n sional geometry can be described (see 
e.g. ICropper et al J 1 1 999t ISuleimanov et al.l [20051) by the mass 
continuity equation 



pV - a. 



(6) 



:7-(pV) = 0, 

dz 

the momentum equation 

dz (z + R^aY 

the energy equation 
dP dV 
az az 
and the ideal-gas law 
pkT 



P = 



(1) 



(2) 



(3) 



(4) 



Here me is the hydrogen mass and y=5/3 is the adiabatic in- 
dex. The cooling rate A due to thermal optically thin radiation is 
given by 



where the mass accretion rate a per unit area is a free parameter 
of the model. 

A simp le analytical solu tion of Eqs. ([T]l-(|4li can be ob- 
tained (see iFrank et al] |2002|) . if one assumes a constant gas 
pressure throughout the PSR and considers the cooling only by 
bremsstrahlung. The distributions of the velocity, temperature, 
and density along the PSR height z in this case are 



V^Vo 



2/5 



T^Tq 



2/5 



P =P0 



3.94 X 10^ a"' fi-^'^ nlm^'^Rf^ cm. 



where the shock height above the surface, 
H, 

and the parameters at the shock position are 

Vo = 1.29 X lO'^m'''^ 
po = 7.75 X 10" 
To = 6.03 X l(f iimR^^ K. 



-2/5 



■RT^^ cms"'. 



am-'I^R]'^ gcm- 



(7) 



(8) 



(9) 



Here yu is the gas mean molecular weight (in units of me) and 
yKe is the number of nucleons per electron in the fully ionized 
plasma (yu - 0.62 and /ie = 1-2 for solar composition). We use 
the dimensionless variables for the WD mass m - M^a/Mq and 
the radius R() = 7?wd/10^cm. The WD radius can be calculated 
from the WD mass-radius relation dNauenberd 1 9721) : 



R9 = 0.78 



[ m ] 1 1.44/ 



1/2 



(10) 



Alternatively, for masses m in the range 0.5-1.2, one can use a 
linear approximation to the mass-radius relation 



R<) = 1.364 -0.807m. 



2.2. Numerical sclieme 



(11) 



A - n^riih^iT), 



(5) 



We now describe the method of numerically solving the PSR 
structure. First of all, the input parameters, the WD mass, and the 
local accretion rate a are defined. Then we use two-loop iteration 
scheme. In the outer loop the amplification factor is changed. 
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while in the inner loop the PSR structure and the emergent spec- 
trum are computed with the fixed distribution of the amplifica- 
tion factor on the relative height z/zo, where zo is the shock posi- 
tion. First, an initial model without Comptonization (A(z) - 0) is 
calculated. Equations O-dlll are solved by the shooting method 
from z = Zo to z = with the following bou ndary conditions at 
z = Zo (see details in lSuleimanov et alll2005h : 



Vo 
po 
Po 
To 



0.25 V2GMwd/(^wd + zo), 
a 

Vo' 
3aVo, 



(12) 



The shock position zo is found iteratively from the additional 
boundary condition: V = at the WD surface. The obtained 
PSR structure along z-coordinate is then considered as a plane- 
parallel atmosphere. The radiative transfer equation in this at- 
mosphere is solved at a grid of 90 column densities E (defined in 
usual way dS - -pdz), distributed equidistantly from E » 10"'' 
g cm"^ to Smax - ^(Rwd)- The Compton amplification factor 
A as a function of the PSR height is given by the ratio of the 
radiation energy loss rate, obtained from the solution of the ra- 
diative transfer equation, to the one given by the thermal cooling 
function: 



dF/dz dF/dl 

1 + A(z) = ; = -p- 



A 



A 



(13) 



where F is the integrated over frequency radiation flux. 

Then we again solve Eqs. ([TJ-lllll with the computed am- 
plification factor Here we assume that the amplification factor 
only depends on the relative position z/zo in the PSR. It is im- 
portant, because in the inner loop of iterations the shock posi- 
tion and the set of geometrical depth z change from iteration to 
iteration. Therefore, for the given shock position zo, the ampli- 
fication factor at height z is computed by interpolating the de- 
pendence A(z/zo). This iteration procedure is repeated until the 
relative change in the shock position is less than 1 per cent (with 
5-6 iterations being sufficient). As a result we obtain the self- 
consistent PSR model with Compton cooling, together with the 
emergent spectrum of radiation. 

2.3. Radiative transfer witli Compton scattering 

The radiative transfer equation for specific intensity I{x, jj.) with 
Compton scattering taken into account in the plane-parallel ap- 
proximation has the form 



dl(x, rj) 



[a(x) + k{xW{x,rf)-S{x,rf)l 



(14) 



where cr{x) and k{x) are electron scattering and true absorption 
opacities at depth S, rj is the cosine of the angle between the nor- 
mal to WD surface and the direction of the radiation propaga- 
tion. The absorption opacity is the sum of the free-free opacities 
of the fully ionized 25 most abundant chemical elements with 
solar composition. The source function is the sum of the thermal 
part and the scattering part: 



S {x, rf) 



k(x) 



cr(x) + k(x) 



B(x) + 



(15) 



x^ djci r r 

~r\ — 77^ ~r \ \ <ij]iR{x,xi,^)I(xur]i), 
o-(x) + k(x) Jo xj J J 



0.1 
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Fig. 2. Dependence of the Compton parameter icomp for the 
PSR on the WD mass. The solid curve giv es the results for the 
analytical model (|7]i of iFrank et alj (|2002|) given by Eq. ( I^TT i. 
while the dotted and dashed curves are for the numerical model 
of Suleimanov et al. (2005) for two accretion rates a - I and 10 
respectively. 



-1 -2 

g s cm 



where B{x) is the Planck function. The Compton scattering 
redistribution function that describes the electron scattering 
(in the isotropic Thomson approximation in the electron rest 
frame, which is ac curate enou gh for electron temperatures below 
100 keV), is ( Arutvurivan & Nikogosvan.1980; .Poutanen. 19941: 
iPoutanen & Svenssonlll996h : 



R{x, xi,tl/) — 



1 



where 



xi = 



hvi 



hv 



(16) 



(17) 



are dimensionless photon energies before and after scattering. 



= 



kT 



(18) 



is the electron temperature in units of the electron rest mass, 

4i - rjriy + Vl ~ ~ '?! '-OS V is cosine of the scattering 

angle, K2 is the modified Bessel function, and 



7* = 



Q 



-y/2xxi(l - if/) 



- x^ + x\ - Ixxii//. 



(19) 



The outer boundary condition is found from the lack of in- 
coming radiation at the PSR surface, and the inner boundary 
condition is taken as the Planck function with the effective tem- 
perature of the WD (« 2 X 10"* K). The formal solution of 
the radiation transfer e quation ([T4| is obtaine d using the short- 
characteristic method (l01son&Kunaszlll987l) . and the full so- 
lution is found with a simple A-iteration method. This method 
gives a fast solution in this case, as the PSR has a rather small 
electron scattering optical depth (t^ < 1). This is also why we do 
not need to consider the nonlinear terms (i.e. induced scattering) 
in the radiation transfer equation. 
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10"^ 10"^ 10"' 1 10 100 

z//?TO E (keV) 

Fig. 3. Comparison of PSR models computed for Mwd = 1.2Mq with a local accretion rate a - 10 g s"' cm"^ with (solid curves) 
and without (dashed curves) Compton scattering taken into account. The temperature and the density structures (panel a), as well as 
the emergent spectra (panel b), are shown. The dotted lines correspond to the analytical solution 



2.4. Results of calculations 

Before presenting the results of our calculations, let us first eval- 
uate the importance of cooling by Compton scattering in the 
PSR. The Compton parameter 



Yc 



omp 



: 40 max(re,Tg) 



(20) 



can be calculated for various models of the PSR. The analytical 
model (|7]) yields 



f 

Jo 



«e dz = 0.05 /z'''^ ^i, 



(21) 



It is clear the Compton parameter can be close to unity only for 
high-mass WDs. It is possibl e to also compute ycomo using nu- 
merical models of the PSRs dSuleimanov et al.ll2005h . The de- 
pendence of icomp on the WD mass is shown in Fig. |2l We see 
that at high accretion rates and high WD masses, Compton scat- 
tering can affect the structure of the PSR and the emergent spec- 
trum. 

The influence of Comptonization on the temperature and the 
density structure of the PSR, as well as the emergent spectra, are 
demonstrated in Fig. [3] Here we take a 1.2 Mq WD with a local 
accretion rate a - 10 g s"' cm"^ and compute the PSR structure 
with and without Comptonization taken into account. The am- 
plification factor 1 + A due to Compton scattering is about 1.5 at 
the top of the PSR (near the shock) and decreases to 1 at the WD 
surface. Therefore, if Comptonization is taken into account, the 
cooling rate becomes higher, the height of PSR smaller, and the 
emergent spectrum slightly harder But the difference between 
the emergent spectra with and without Comptonization is, how- 
ever, rather small. 

We computed the grid of PSR models for WD masses rang- 
ing from 0.3 to 1.3 Mq with step 0.02 Mq (51 models in total). 
The total accretion rate was assumed to be M = 10'^ g s"' with 
the fixed part / = 10"^ of the WD surface occupied by PSR. This 
gives different local accretion rates a for various WD masses, 
varying from ~ 0.5 to ~ 5 g s"' cm"^. The emergent spectra of 
PSR with different local acc retion rates (in the inte rval 0.1-10 
g s ' cm"^) are very similar dSuleimanov et al.ll2005h : therefore, 
we have only one free parameter, the WD mass, in our grid of 
model spectra. This means that the mass can be found by fitting 



the observed X-ray spectrum. Of course, the WD mass deter- 
mination using spectral fitting is model dependent, because of 
using of the specific theoretical WD mass-radius relation and of 
the framework the theoretical spectra are calculated. However, 
the WD mass determination depends very little on the particular 
choice of M and /. 

The grid of the IP spectral models has been incorporated into 
the XSPEC software package under the name IPCOMP and is 
publicly available. 



3. X-ray spectrum of V709 Cas 

V7 09 Cas is a typical IP emitting in the hard X-rays (see 
e.g. ^Suleimanov et al.ll2005b . The pre sent dataset was ob tained 
using the RXTEfPCh (3-20 keV) dJahoda et all 1 1996) and 
/Afr£GiML/IBIS/I SGRI (20-120 keV) (Ubertini etaOlzOOSl; 
iLebrun et al.ll2003l) instruments with the total effective exposure 
times of 47 ks and 4.5 Ms, respectively. Our spectral model con- 
sists of the PSR model IPCOMP, an iron line modelled as a 
Gaussian, and the partial absorption model PCFABS described 
by the covering factor Cp and the column de nsity A^h. The spec- 
tral analysis is performed with the XSPEC dArnaudll 19961) ver- 
sion 12.4. 

The best fit with ;t'^/d.o.f. = 31/30 is obtained for Mwd = 
O.9l![J [J^M0, a partial covering fraction Cp — 0.41 + 0.05 and 
A^H = (53 + 10) X 10^^ cm"^ (see Fig.|4|l. The iron line energy 
is 6.5 + 0.2 keV with an equivalent width of 508 eV. The un- 
absorbed flux between 0.5 and 150 keV is F ^ 2.4 x 10"'° 
erg cm"^ s"', which translates into a bolometric luminosity of 
L ^ 1.52 X 10^^ erg s~', for the source distance of 230 pc 
dBonnet-Bidaud et alj|2()0lh . 

We can estimate the mass accretion rate using the simple re- 
lation M = LR^JiGM^i). We get M = 0.8 x lO"" g s"' for 
the WD radius of /?wd - 0.61 X 10^ cm, which was computed 
using Eq. dTOl i. To investigate the role of Compton scattering, 
we also implemented into XSPEC the grids of the PSR model 
without Compton scattering. The derived parameters are simi- 
lar to those obtained with the Compton scattering model with 
a comparable reduced The WD mass is found to be slightly 
lower Mwd = O.88;^g gjM0. This seems counter-intuitive, because 
Comptonization increases the flux at high energies, and there- 
fore would need lower temperature (and therefore lower WD 
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Fig. 4. The unfolded spectrum of V709 Cas and the best-fit spec- 
tral model of the PSR. The crosses represent the data from 
RXTEfPCA (3-20 keV) and INTEGRALflSGRl (20-100 keV) 
instruments. The total spectrum of the PSR model is shown 
by the solid (with Compton scattering) and the dashed (with- 
out Compton scattering) curves. The dotted curve represents the 
Gaussian iron line. 



mass) to radiate the same energy. However, the spectral slope 
in the 4-20 keV energy range is different for the two models, 
so the absorber parameters also change to Cp - 0.44 + 0.05 
and A^H = (63 + 10) x 10^^ cm^^ Thus for V709 Cas, both 
mass estimations are consistent with each other within the er- 
rors, and Compton scattering does not have much effect, which 
is expected for such a WD mass and a rather low accretion rate. 

For illustrative purpose, we also obtained a fit of the 
observed spectrum with a single temperature bremsstrahlung 
model. We obtained T = 29.6 + 2.5 keV, Cp = 0.31 + 0.04, 
A^H - (58 + 15) X 10^^ cm"^, and the iron Une energy 6.5 + 0.2 
keV with the reduced = 1.062. Using a linear approximation 
to the mass-radius relation (fTTT i. we can estimate the WD mass: 



T'keV 



23.6 + 0.59rke 



(22) 



which gives M„d = 0.72 + 0.035Mo. This result illustrates the 
error from the use of a simple one-temperature bremsstrahlung 
model for the WD mass estimation. 



4. Conclusions 

We have studied the influence of Compton scattering on the 
structure of the post-shock region and the broad-band X-ray 
spectra of intermediate polars. Compton scattering leads to 
slightly harder emergent spectra and smaller PSR height. The ef- 
fect can be significant for the luminous intermediate polars with 
high-mass WDs. We have also constructed a grid of the emergent 
spectra in a wide range of WD masses and incorporated into the 
XSPEC package. We used this grid to fit the broad-band X-ray 
spectrum of V709 Cas observed by RXTE and INTEGRAL obser- 
vatories and determined the WD mass in this object. Accounting 
for Compton scattering does not change the obtained mass sig- 
nificantly. 

In this work we have not considered the cooling by cyclotron 
radiation. However, for WDs with stronger magnetic field (po- 
lars), Comptonization of the cyclotron photons can significantly 
affect the structure of the PSR and the emergent spectra. We plan 
to perform such a study in the future. 
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